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In this paper we present a review of bosonic renormalization eflects on electronic carriers observed 
from angle-resolved photoemission spectra in the cuprates. We specifically discuss the viewpoint 
that these renormalizations represent coupling of the electrons to the lattice, and review how the 
wide range of materials dependence, such as the number of Cu02 layers, and the doping dependence 
can be straightforwardly understood as arising due to novel electron-phonon coupling. 



I. INTRODUCTION 

The discovery of a "kink" in the nodal ((0,0) 
- (7r,7r)) dispersion of the high-Tc cuprates, and 
band renormalizations, in the form of a peak-dip- 
hump structure in the anti-nodal (0,7r) - (7r,7r) 
dispersioni^>^^>^>&i>a>a^ii^i^^i^^ii, have at- 
tracted considerable attention in recent years. These 
band renormalizations are interpreted as due to electron- 
boson coupling, and it is believed that understanding 
their origin will provide information about the underly- 
ing pairing mechanism in these materials. There is still 
considerable debate as to the identity of the responsible 
bosonic modei^. 

One proposal is coupling to an electronic mode associ- 
ated with the collective mode found in neutron scattering 
near momentum transfers (7r,7r), the so called magnetic 
resonance mode^i^ i^^i^^ . Additionally, renormalizations 
could also be due to coupling of electrons to damped 
magnons, which have less well-defined momentum struc- 
ture. However, the strength of the coupling of magnons 
to electrons is under debatei^i^. For example, quan- 
titative comparisons of ARPES and neutron measure- 
ments on YBa2Cu306.6 (YBCO) have been made, and 
the overall strength of the coupling inferred from the data 
was indicated to be of sufficient strength to give rise to 
superconductivitjiii. We also remark that a quantitative 
comparison between the neutron scattering and ARPES 
measurements reported in Ref. [13 can be complicated 
by the polar surface of cleaved YBCO resulting in sur- 
face reconstruction. This can produce significant differ- 
ences between the bulk and surface layers of this mate- 
rial. Moreover, a systematic examination of the effects 
of doping, material class, and temperature have not been 
thoroughly explored. In fact, given that the spin contin- 
uum arises from the Cu02 plane, one might expect the 
coupling to electrons to be relatively material-class inde- 
pendent. On the other hand, it is well-known that a neu- 
tron resonance displays a material dependence, appear- 
ing at larger energies for larger Tc materials including 



both the single and multi-layer cuprates. This opens the 
possibility of linking the neutron resonance with ARPES 
renormalizations via a material-dependent study. 

An alternative proposal is coupling to a spectrum 
of oxygen vibrational phonon modes^i^ i^^i^^ , specifically, 
the c-axis out-of-phase bond-buckling oxygen vibration 
or Big mode ($7 35 — 45 meV) and the in-plane bond- 
stretching oxygen mode (fJ ~ 70 — 80 meV). This pro- 
posal has been able to account for many experimental 
observations including the anisotropy of the observed 
renormalizations^!, fine structure in the form of sub- 
kinks observed in the temperature dependence of the self- 
energy in Bi2Sr2CaCu208+^ (Bi-2212), and doping de- 
pendent changes in self-energ y^^i^^ . This interpretation 
is further supported by recent ARPES experiments that 
have measured an ^^O ^^O isotope shift in the nodal 
kink position^. Finally, analogous features have been 
observed in tunneling spectra which have also been in- 
terpreted as coupling to a bosonic mod o^^'^^i^^'^^ , and 
possibly the same mode responsible for the ARPES ob- 
served kink. Ref. [2^ reported an ^^O isotope shift of the 
feature in a Bi-2212 sample indicating a phonon origin. 

In this work, we present a review of how these two 
scenarios can be differentiated by studying the material, 
doping and temperature dependence of the band renor- 
malizations. Our focus is on how c-axis phonons provide 
a material dependence to the ARPES kinks due to the 
phonon's sensitivity to local symmetry and the environ- 
ment surrounding the Cu02 plane. We consider doping 
dependent changes to the renormalization in Bi-2212 as 
well as the dependence of the renormalization within the 
Bi and Tl families as the number of Cu02 layers is var- 
ied. We will also discuss some recent ARPES results 
on the n = 4 layer system Ba2Ca3Cu408F2 (F0234)^. 
In this system, the inner and outer layers occupy differ- 
ent crystal environments resulting in differing Madelung 
energies associated with each plane in the undoped com- 
pound. This difference drives inequivalent dopings be- 
tween the two sets of layers, with one set n-type and 
the other p-type. The inequivalent doping in each plane 
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generates further symmetry breaking in the layers and 
the el-ph coupUng is expected to differ. Indeed, Ref. [1^ 
observes stronger kink features in the plane associated 
with the outer (n-type) layer of the material. Here, we 
will discuss these observations in the context of the el-ph 
coupling scenario. 

The organization of this paper is follows. In section II 
we discuss the role of symmetry breaking in producing 
el-ph coupling to c-axis phonons and how such coupling 
is expected to vary with the crystal structure of the high- 
Tc cupratcs. In section III we present ARPES data for 
various multi-layer cuprates in order to examine how the 
now the band renormalizations vary with the number of 
Cu02 layers, n. Energy distribution curves in the anti- 
nodal region of the Bi and Tl-families for n = 1 — 3 and 
MDC dispersions for the p-bands of the F-family with 
71 = 4 — 5 are presented. Of particular interest are the 
results for the single layer Tl cuprate Tl-2201. Here, we 
show that Tl-2201 does not resolve the typical peak-dip- 
hump structures in the antinodal region despite the fact 
that the spin resonance mode exists in this system^. In 
section IV the doping dependence of the nodal and antin- 
odal dispersions for Bi-2212 are presented. The renormal- 
izations in each region behave differently as the samples 
are overdoped, pointing to presence of multiple bosonic 
modes. In section V we present a theoretical basis for 
understanding the self-doping phenomena in the parent 
compound F0234 and discuss what implications this pro- 
cess has on the coupling to c-axis phonons. Finally, in 
section VI, we conclude with a brief summary and some 
additional remarks. 



II. ELECTRON-PHONON COUPLING IN 
MULTI-LAYER CUPRATES 



In this section we discuss how c-axis phonons can 
be sensitive to the material environment off the Cu02 
planes. Electron-phonon (el-ph) coupling to the in-plane 
bond-stretching modes are of a deformation type and the 
strength of the couphng depends on the Cu-0 bond dis- 
tance. Since this distance is relatively constant in all 
cuprates, coupling to these modes is relatively materials 
independent. This is not true for coupling to the c-axis 
bond buckling modes. Since our goal is to explore the ma- 
terials dependence of the band renormalizations, we will 
focus our attention to the c-axis modes. El-ph coupling 
to c-axis phonons can arise to first order from modula- 
tions in the on-site energy due to atomic site oscillation 
through local crystal fields. For the Cu02 planes, such 
crystal fields arise due to asymmetries in the environment 
surrounding the plane2ii2£ii^ and the degree of symmetry 
breaking varies from material to material with the chem- 
ical environment (number of layers and doping from the 
ideal stoichiometric compound). 

In an ideally undoped single layer cuprate the structure 
above and below the planes is identical and little sym- 
metry breaking occurs across the plane. The coupling to 
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FIG. 1: The local crystal field strength at the planar oxygen 
site of the outermost Cu02 plane of the Hg-family of cuprates. 
All results have been normalized by the maximum field which 
occurs for the n = 3 layer system. 



the planar c-axis modes is therefore expected to be weak 
in the parent compound. In the undoped multi-layer sys- 
tems symmetry breaking across the outer planes can be 
large leading to a sizeable coupling in these layers. A sec- 
ond pathway for symmetry breaking occurs as the parent 
compounds are doped by introducing substitutional or in- 
terstitial dopant atoms in the charge reservoir area off the 
Cu02 planes. These dopants donate charge to the Cu02 
plane, and cast off-plane electric fields that cannot be 
effectively screened by the in-plane carriers^. As a result 
they form charge impurities that further break symmetry 
across the Cu02 planes. Through this mechanism, cou- 
pling to c-axis phonons can occur in single layer systems 
where they are normally forbidden by symmetry. 

To quantify the effect of the crystal environment in the 
parent systems, we calculate the local crystal field at the 
outermost Cu02 layer for the HgBa2Ca„Cu„02(„+i)+5 
(n = 1-6) family of cuprates. To do this, we use ex- 
perimental structural dataH and assume an ionic point 
charge model with formal valences assigned to each atom. 
The Ewald summation method^ is then used to perform 
the electrostatic sums for the electric field strength. The 
results are shown in Fig. [1] 

In the single layer compound the local field is zero since 
the Cu02 layer coincides with the plane of mirror sym- 
metry and we have not considered dopings away from 
pure stoichiometry. In the double- and triple-layer com- 
pounds the value of the field rises taking a maximum 
for n = 3. After this, as the number of layers further in- 
creases, there is a general trend of decreasing field values. 
(We note here that the structural data for the n = 5 and 
n = 6 compounds had a large degree of error, presum- 
ably from the difficulty in sample growth.) This decrease 
is due to a reduction in the degree of symmetry break- 
ing across the outer plane once the number of layers is 
increased beyond n = 3^. 

In the limit of the infinite layer system CaCu02 the en- 



vironment is symmetric around every plane. In this case 
the electrostatic contribution to the coupling is identi- 
cally zero however steric forces introduce buckling to the 
plane which creates a new pathway for el-ph coupling 
of a deformation type^^. However, in terms of material 
dependence, such steric forces are present in all Cu02 
systems and therefore do not contribute to differences 
between materials. 

The strength of the coupling to c-axis modes, such as 
the Big mode, scales as A oc E'^. Furthermore, due to 
the oxygen charge transfer form factos, the Bi^ mode 
couples most strongly to antinodal electrons^!. Therefore 
the behavior of the antinodal renormalization is expected 
to have a dramatic materials dependence if it is due to 
the Big phonon. This layer dependence is expected to be 
different for coupling to the spin-resonance mode. There- 
fore, by examining the change in coupling strength as one 
moves up in the number of layers provides a pathway for 
distinguishing between these two scenarios. 



III. LAYER DEPENDENCE 

In the previous section we discussed how the material 
dependence of the renormalizations is expected to arise 
in the various families of cuprates and how this coupling 
is expected to differ in the phonon and spin resonance 
proposals. We now wish to review the available ARPES 
data in light of the theoretical considerations of the pre- 
vious section. Here, our focus is on the observed changes 
in the renormalizations as the number of layers within 
the Bi- and Tl-families. 

Single crystals of nearly optimally doped 
TlaBasCaCuzOg (Tl-2212), TlBazCaaCugOg (Tl- 
1223) and slightly overdoped TlzBazCuOe (Tl-2201) 
were grown using the flux method. As-grown T12212 
(Tc = 107 K) and Tl-1223 (T^ 123 K) crystals were 
chosen for the ARPES measurements. Tl-2201 crystals 
used in our measurement were prepared by annealing 
the as-grown crystal (Tc ~ 30 K) under a nitrogen flow 
at a temperature of 500°C, yielding a T^ of 80 K. The 
data were collected using a Scienta R4000 photoelectron 
spectrometer. Measurements were performed at the 
SSRL beam line 5-4 using 28 eV photons and at the 
Advanced Light Source beam line 10.0.1 using 50 eV 
photons. The energy resolution was set at 15-20 meV 
for the Tl data presented in this work. Samples were 
cleaved and measured in ultrahigh vacuum (< 4 x 10~^^ 
Torr.) to maintain a clean surface. Detailed ARPES 
results on these compounds have been reported in Rcf. 

M 

Although the dispersion kink along the nodal direc- 
tion has been found universally in high-Tc cuprates^, the 
momentum dependence of this renormalization feature, 
when moving away from the nodal direction, exhibits a 
material depcndence^i^i. It has been confirmed recently 
that there is dependence on the number of Cu02 planes 
in the unit cell in the Bi-family and, most recently, in the 
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FIG. 2: Representative EDCs near the antinodal region of the 
Bi- and Tl-families of cuprates, including single layer (Bi-2201 
and Tl-2201), bi- layer (Bi-2212 and Tl-2212), and tri-layer 
(Bi-2223 and Tl-1223) compounds. The high background in 
the data of Tl-1223 is probably due to the absence of a natural 
cleaving plane in the crystal structure. Nevertheless, a peak- 
dip-hump structure in the spectrum can still be discerned. 
The red dashed line is a guide-to-the-eye to make the "hump" 
more discernible. 



Tl-family of cuprates^. In the multi-layer compounds, 
the "kink" becomes more dramatic and eventually breaks 
the band dispersion into two branches: one branch with 
a sharp peak and another branch with a broader hump 
structure^iiiiSi^i^l. The two branches asymptotically ap- 
proach one another at a characteristic energy scale of 
70 meV and coincide with the dominant energy scale of 
the kink along the nodal dispersion for nearly optimally- 
doped cuprates. This separation of the band dispersion 
becomes most prominent near the antinodal region and 
results in the famous peak-dip-hump structure^i in the 
energy distribution curves (EDCs), as shown in Fig. ^ 

The momentum dependence of the kink is quite dif- 
ferent in the single layer compound, where the disper- 
sion kink becomes less prominent moving away from the 
node. In addition, the band dispersion retains a single 
branch with no separation observed, unlike the case of 
the multi-layer compounds^i^^i^. As a result, no appar- 
ent peak-dip-hump structure can be seen in the EDCs 
near the antinodal region for the single layer compounds 
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Fig. sa. 

In summary, the layer dependent renormalization near 
the antinodal region is due most hkely to electrons cou- 
pled to a sharp bosonic mode, whose origin is strictly 
constrained by the number of layers in the material. This 
mode is either absent, or has a negligible coupling to the 
electrons, in single layer compounds, but exhibits promi- 
nent coupling in the multi-layer compounds. The spin 
resonance mode does exist in some single layer cuprates^ 
(notably Tl-2201 by not La2Cu04+5). Therefore, one 
can conclude that the spin resonance mode is an unlikely 
candidate for the mode responsible for the renormaliza- 
tions in the antinodal region. On the other hand, cou- 
pling to c-axis phonons can exhibit a very different cou- 
pling in single- and multi-layer compounds. As we have 
discussed, the Big phonon couples strongly to the elec- 
trons in multi-layer compounds and weakly to electrons 
in single layer compounds. This mode can also reproduce 
the observed anisotropic momentum dependence of the 
renormalization in bi-layer Bi-22132i. We also note that 
the form for the Big coupling is attractive in the d-wave 
pairing channel^, which could be one factor enhancing 
Tc in the multi-layer systems. 

We now consider another case, that of the F-family 
of cuprates with n Cu02 layers, n = 3 — 5. The single 
crystalline samples were grown by the flux method under 
high pressure^. ARPES measurements on the F-family 
were performed at beamline 10.0.1 of the Advance Light 
Source (ALS) at Lawrence Berkeley National Laboratory. 
The measurement pressure was kept < 4x 10^^ Torr at all 
time and data were recorded by Scienta R4000 Analyzers 
at 15K sample temperature. The total convolved energy 
and angle resolution were 16 meV and 0.2° respectively 
for photoelectrons generated by 55 eV photons. 

In Fig. [3] we present MDC derived dispersions for p- 
type bands of the three (F0223), four (F0234) and five 
(F0245) layer F-based cuprates. In all three cases the 
dispersions show clear kinks, but at a larger energy scale 
in the four- and five-layer materials. In a d-wave super- 
conductor coupled to an Einstein mode, the energy scale 
of the kink occurs at 51 -I- Aq where il is the energy of 
the mode and Aq is the maximum value of the supercon- 
ducting gap. This shift in energy scale is a reflection of 
the changes in the superconducting gap size as n is varied 
from 3 to 5. In order to quantify the strength of the kink, 
slopes are extracted from the dispersion above and below 
the kink position, de/dk\y and de/dk\^. An estimate for 
the relative coupling stengths A' is then given by; 



This procedure produces A' = 0.89, 0.75 and 0.49 for 
n = 3, 4 and 5 respectively. This trend is easily under- 
stood in the phonon mode scenario where the dominant 
mode in the superconducting state is the Big mode for 
which the coupling strength is proportional to the local 
crystal field. The observed decrease in coupling strength 
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FIG. 3: MDC derived dispersions along the nodal direction 
(0,0) - (7r,7r) of thep-type band in the 3-layer (F0223), 4-layer 
(F0234) and 5-layer (F0245) F-family of cuprates. 

can be easily understood from the expected local field 
strengths in the outer most layers (recall Fig. [T]), which 
decreases for n > 3. 



IV. ENERGY SCALES AND DOPING 
DEPENDENCE 

In the el-ph coupling picture the carriers couple to a 
spectrum of bosonic modes and we have already seen how 
the c-axis modes can produce a materials dependence of 
the renormalizations. It is important to note that the 
coupling to each of these modes is highly anisotropic 
with the in-plane bond-stretching (breathing) mode dom- 
inating in the nodal region and the out-of-plane bond- 
buckling (Big) mode dominating in the anti-nodal region. 
Since these modes have different frequencies one would 
naturally expect the multiple energy scales to manifest 
in the experimental data. Indeed, evidence for multiple 
energy scales has been found both in the temperature 
dependenceii2ii^ as well as the doping dependence of 
Bi-22132ii^i^. In this section we revisit the doping de- 
pendence of the nodal and antinodal renormalizations, 
highlighting the different behavior in each region of the 
Brillouin zone and discuss how this dichotomy further 
supports the el-ph scenario. 

High quality single crystals of optimally doped 
Bi2Sr2Cao.92Yo.o8Cu208+5 (Bi2212 OP, T, = 96 K) were 
grown by the floating zone method. The overdopcd crys- 
tals with Tc = 88 K were prepared by post annealing 
the optimally doped Bi2212 crystal under oxygen flow 
at a temperature of 400° C. The overdoped sample with 
Tc = 65K is a derivative of the Bi2212 family with lead 
doped into the crystal to achieve such an overdoped con- 
figuration. The data were collected by using He I light 
(21.2 cV) from a monochromatcd and modified Gamma- 
data HE Lamp with a Scicnta-2002 analyzer and in SSRL 
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FIG. 4: The doping dependence of the antinodal spectrum of 
Bi2212 taken in the superconducting state (lOK). Shown in 
the upper row are the false color plots of the spectra taken 
along the indicated cut direction (inset). The black dots are 
the peak and hump positions of the bonding band seen in the 
EDCs. Shown in the lower row are the EDCs along the dashed 
line indicated in the false color plots. The symbols "AB" and 
"BE" represent the antibonding and binding bands while the 
numbers are the energy position of the dip of the EDC. 
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FIG. 5: The doping dependence of the nodal spectrum at a 
temperature well below Tc (10 K). Shown in the upper row 
are the false color plots of the spectra taken along the cut 
as indicated in the inset. The black curves are band disper- 
sion obtained by fitting MDCs to Lorentizan functions. The 
apparent kink position in the dispersion are marked by the 
yellow shaded area, which appears to be approximately the 
same for all three dopings. The dashed lines serve as a guide- 
to-the-eye for visualizing the apparent kink in the dispersion. 
Shown in the lower row are the real part of the self-energy 
extracted from subtracting the band dispersion from a linear 
bare band. The arrows indicate the positions of fine structure 
in the extracted ReE. 



beamline 5-4 using 19 eV photons with a Scicnta-200 an- 
alyzer. The energy resolution is '^lO meV and angular 
resolution ~0.35°. The samples were cleaved and mea- 
sured in ultra high vacuum (< 4 x 10~^^ Torr.) to main- 
tain a clean surface. 

In Fig. [4] ARPES data taken along a cut in the antin- 
odal region illustrate this effect. The upper panels show 
the measured spectral function along the cut while the 
lower panels show ^(k, uj) at a fixed k-point as indicated 
by the dashed lines. 

Near (0,7r), the energy of the dip feature is the best 
measure of the energy scale of the mode responsible for 
the renormalization^iiii. For the optimally doped sample 
(OP96K), shown in the first column of Fig. 21 the dip po- 
sition is clearly located at 70 meV. This can be seen 
in both the false color plot and the EDC cut. For moder- 
ate overdoping (OD88K), the energy of the dip is lowered 
to ~ 58 meV while for heavily overdoped (OD65K) the 
energy is lowered further to ~ 32 meV. In both over- 
doped cases, contributions from the bonding- (BB) and 
anti-bonding (AB) bands contribute to the quasiparticle 
peak at the Fermi level. In the OD65K case, the con- 
tribution from the AB makes an exact determination of 
the dip position difficult and the estimate of ~ 32 meV 
should be considered a lower bound. 

Turning now to the nodal region, we find qualitatively 
different behavior. Fig. [5] presents A(k, w) along the 
nodal cut ((0,0) - (7r,7r)) for the same three samples. The 



highlighted region indicates the approximate position of 
the kink. In the nodal region the overall bandwidth is 
much larger than the energy of the bosonic modes so 
the dramatic band breakup does not occur— and the 
renormalization manifests as a kink in the dispersion. 
In this case, the energy scale of the kink is most eas- 
ily determined from the structure of the real part of the 
self-energy ReE. The MDC-derived estimate for ReE, 
obtained from subtracting the MDC-derived dispersion 
from an assumed linear band, is also in the lower panels 
of Fig. El 

The nominal doping dependence^ of the energy scales 
in the nodal and antinodal region of Bi2212 are sum- 
marized in Fig. in^. For reference, the superconducting 
gap Aq is shown also, which is determined from the peak 
positions of the Fermi function divided spectrum at the 
Fermi level. While the characteristic energy in the antin- 
odal region (dip energy) follows the decrease in the su- 
perconducting gap, the characteristic energy in the nodal 
region remains more or less constant (~ 70 meV). The 
difference in the doping dependence of the two energy 
scales lends further support to the existence of coupling 
to multiple modes. If a single mode were responsible for 
the renormalization throughout the zone one would ex- 
pect the doping dependence to follow the same trend in 
the nodal and antinodal regions. 

The energy of the dominant mode can be obtained 
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FIG. 6; A summary of some energy scales relevant to the 
renormalized band dispersions, (a) The apparent kink po- 
sition in the nodal dispersion, dip energy at the antinodal 
region and the superconducting gap are summarized for the 
three doping levels shown in Fig |4] (b) the mode energy 
obtained by subtracting the superconducting gap from the 
characteristic energies of the renormalization effect. The red 
arrow is to remind the reader that the shown quantity at the 
antinodal region of the OD65K sample is a lower bound for 
the actual value. 



by subtracting the magnitude of the superconducting gap 
from the observed energy scale, expected to be + A(j^. 
The results of this procedure are shown in Fig. [SJs. The 
energy of the dominant mode in the antinodal region, 
within the error bars of the data, is independent of dop- 
ing. The behavior in the nodal region is different; the 
energy of the dominant mode changes with doping. At 
optimal doping the energy of the dominant mode is ^ 35 
mcV but in the overdoped samples the energy is larger 
^ 60 meV. This result is consistent with the picture of 
coupling to multiple modes outlined in Ref. [2l|. We also 
note that in OD88K a secondary feature can be observed 
in ReS at precisely the same energy as the dip energy of 
the antinodal region. Similar fine structure was reported 
earlier in Ref. |22| . The presence of this sub-feature in the 
UD88K data as the sample is progressively overdoped, 
along with the 35 meV scale in the nodal data at op- 
timal doping, is evidence of a trade off between a cou- 
pling dominated by the Bi^ mode and one dominated by 
the bond stretching mode. We further note that Ref. |8| 
reached similar conclusions but assigned the anti-nodal 
renormalizations to the spin resonance mode. We be- 
lieve that the multi-layer data of the previous section, 
especially the single-layer Tl data which shows no renor- 
malization in the anti-nodal region, directly refutes this 
conclusion and favors the el-ph scenario. 



V. SELF-DOPING 

We now turn our attention to the identity of the car- 
rier types in the inner and outer planes in F0234. Ex- 
perimentally, the parent compound of F0234 is known 
to self-dope with the inner and outer planes having 
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FIG. 7: A summary of the electrostatic calculations for F0234. 
(a) The Madelung potential difference between the Cu (red) 
and O (black) sites of the inner and outer planes as a func- 
tion of the oxygen buckling of the outer layer. The shaded 
region indicates the region where the inner and outer planes 
are identified as p- and n-type respectively. Buckling of the 
planar O away from the mirror plane of the crystal is defined 
to be positive, (b, inset) The filling of the inner (red) and 
outer (black) planes as a function of the site energy difference 
(buckling) between inner and outer layers. Fillings are ob- 
tained with a simple tight- binding model (see text), (b) The 
electric field at the oxygen site of the inner (black) and outer 
(red) layers as a function of the buckling distance. 



Fermi surfaces of different carrier type. A recent ARPES 
studjs2£ found that the p-typc bands are bilaycr split 
along the nodal direction. Since the inner planes are ex- 
pected to have a stronger inter-planar coupling than the 
outer planes, the observation of bilayer splitting provides 
strong evidence that the inner layers are p-type while the 
outer layers are n-type. Furthermore, kinks in the nodal 
dispersion of both sets of planes but the strength of the 
coupling in the n-type layer by a factor of two. 

The doping of individual layers will be driven by 
Madelung energy differences between the inner and outer 
planes. In order to determine the Madelung potential 
<i> of each site we again employ the Ewald summation 
technique using the structural data from Ref. [4^ and 
assigning formal valence charges to each atom. In princi- 
ple, both the inner and outer Cu sites, as well as the 



7 



inner and outer O sites, can have different Madclung 
energies. The in-plane difference A$" = $q — <i?^„, 
where a is a plane index, determines the charge transfer 
energy^ A" for each plane. The difference between lay- 
ers = <I>°^ — <i>*^, where i = Cu or O, controls the rel- 
ative site energies of the two planes. However, our Ewald 
calculations show that A*p ~ A°p and A4>^„ ~ A$^. 
We therefore neglect the differences and take the charge 
transfer energy to be the same in each layer but shift the 
outer plane's site energies by Ae""- = A$/e(cx)) where 
e(c») is the dielectric constant which taken to be 3.5. 

To understand the self-doping phenomena we consider 
a simple tight-binding model for the Cu02 planes. For 
simplicity, the form for the bandstructure is taken from 
a 5-parameter tightbinding fit to the low-energy disper- 
sion of Bi2212iSi as a representative singleband model. 
We also introduce the usual inter-planar coupling terro^ 
ej^(k) = tperp{cos{kxa) — cos{kya))'^ /A with t± = 50 meV. 
The resulting model Hamiltonian is 

+ E<„^„'>^^(k)[dl_^k^„,.,k + ^-C-] (2) 

where a = 1 — 4 is the plane index, eQ.(k) = e(k) -|- 
A$^/e(oo) for the outer planes and ea(k.) = e(k) for the 
inner planes and < ... > is a sum over neighboring planes. 
(The in-plane charge transfer energy has been absorbed 
into the definition of e(k).) Finally, n is the chemical 
potential which is adjusted to maintain the total filling 
of the parent compound, Ea a- ^o-,a = 4. The resulting 
model is then diagonalized for a given A$ in order to 
obtain the relative filling of the four planes. 

In previous LDA treatments of this material^, it was 
found that the identification of the carrier types in each 
layer is dependent on the degree of Cu-0 buckling oc- 
curring in the outer layers. We therefore calculate the 
Madelung potential difference between the inner and 
outer layers as a function of the degree of buckling. The 
results are shown in Fig. [7^. For flat or dimpled {5 < 0) 
planes, A<i> = — *&ip is positive indicating a larger 
site energy in the outer layer. Charge is therefore ex- 
pected to flow from the outer planes (OP) to the inner 
planes (IP) and the OP/IP are identified as p- and n- 
type respectively. As the plane is buckled outward A<I> 
is suppressed and becomes negative for (5 ~ 0.04 A at 
which point the identity of the planes is reversed. The 
filling of the IP and OP, obtained from diagonalizing 
is shown in the inset of Fig. [7}3. 

Once the fillings have been obtained wc uniformly 
transfer charge between the planes and recalculate the 
crystal field at the oxygen sites. The resulting fields are 
plotted in Fig. [71d. For a dimpled plane the electric field 
at the IP and OP roughly follow a linear dependence on 
the buckling distance with the IP having only a slightly 
weaker field strength. In this case, the kink strength 
measured by ARPES in the two bands would be compa- 
rable. For a buckled plane, once the sign of the Madclung 



potential difference changes and the identification of the 
planes reverses, deviations in the field strength become 
more pronounced. Since the coupling enters as the 
n-type outer plane is expected to have a stronger kink 
feature, consistent with photoemission experiments. 

The carrier concentration of the two sets of planes, de- 
termined from the Luttinger fraction, were reported in 
Rcf. \M with dopings of 0.60 ± 0.04 and 0.4 ± 0.03 in the 
p- and n-type bands. This observation, coupled with the 
different kink strength in each of the bands corresponds 
to (5 = 0.1 A in our calculations. For this degree of buck- 
ling, the ratio of the crystal fields is Eop/Eip = 1.55 and 
the ratio of the kink strength is Xopl\p = 2.4. (The 
absolute value of A depends on the form of the coupling 
constant (^(k, q)2ii^.) This ratio is slightly larger than 
the ratio of ~ 2 reported in Rcf. however the agree- 
ment is good considering the simplicity of the model. It 
is also clear that the incquivalcnt dopings can be under- 
stood from electrostatic configurations. Furthermore, the 
asymmetry in the coupling in each layer falls naturally 
into the c-axis phonon scenario with the strength of the 
coupling being further driven by the self-doping process. 



VI. CONCLUSIONS 

In this work we have presented aspects of the mate- 
rial and doping dependence of the dispersion renormal- 
izations in the nodal and antinodal regions of various 
single- and multi-layer cuprates. We have found that the 
strength of the nodal kink has a strong material depen- 
dence and varies with the number of layers present in the 
material. In general, the kink strength mirrors Tc, taking 
on a maximal value in the n = 3 compounds. The issue 
can be complicated further in the multi-layer cuprates, 
where Madelung potential differences can lead to inequiv- 
alant dopings in the various layers. This can lead to 
further symmetry breaking across the Cu02 planes, and 
results in different kink strength in the different layers 
within the same material. Using a simple tight-binding 
model and electrostatic calculations we have developed a 
picture of this phenomena in self-doped F0234, which is 
consistent with recent ARPES studies. 

The renormalization in the antinodal region also shows 
a marked dependence on the number of layers present 
in the material and is unresolved in the single layer 
cuprates. This result is difficult to reconcile for coupling 
to the spin resonance mode, which is expected not to vary 
with the number of layers, but is naturally explained by 
coupling to the Bi^ phonon when one considers the crys- 
tal structure of these materials. 

Further evidence for multiple phonon modes was found 
in the doping dependence of Bi-2212. Here, the features 
in the nodal and antinodal regions exhibit different be- 
havior. Once gap referenced, the energy scale in the 
nodal region changes from ~ 35 — 40 mcV to ~ 70 — 80 
mcV as the sample is overdoped. This change of energy 
scales cannot be explained by coupling to a single mode 
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and therefore rules out the spin resonance mode, at least 
as the sole player. In the phonon scenario this signifies 
a trade-off between dominant coupling to the Bi^ mode 
near optimal doping and a dominant coupling to the bond 
stretching mode in the overdoped samples. The change 
of relative coupling is due to increased screening of the 
Big mode as the carrier concentration is increased. 

Both the doping and materials dependence presented 
here provides compelling evidence that a spectrum of 
phonon modes are responsible for both the nodal and 
antionodal low-energy renormalizations observed in the 
cuprates. 
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